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Abstract 

We investigate the conditions for a QCD axion to coexist with stabilised 
moduli in string compactifications. We show how the simplest approaches 
to moduli stabilisation give unacceptably large masses to the axions. We 
observe that solving the F-term equations is insufficient for realistic moduli 
stabilisation and give a no-go theorem on supersymmetric moduli stabil- 
isation with unfixed axions applicable to all string compactifications and 
relevant to much current work. We demonstrate how nonsupersymmetric 
moduli stabilisation with unfixed axions can be realised. We finally out- 
line how to stabilise the moduli such that f a is within the allowed window 
10 9 GeV < f a < 10 12 GeV, with /„ ~ ^M SU syM p . 
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A Constraining Corrections to the Kahler Potential 

1 Introduction 

String theory remains certainly the most popular and on all appearances the 
best candidate for an ultraviolet completion of the Standard Model that will 
unify gauge and gravitational interactions in a consistent quantum theory. Since 
the discovery of Calabi-Yau compactifications [1] , it has been hoped that the mea- 
sured properties of the Standard Model may be understood through the details 
of a string compactification. 

This has traditionally been approached top-down, starting with a fully speci- 
fied and globally consistent closed string compactification. The defined technical 
problem is to construct a compactification whose low energy gauge group and 
matter content is that of the Standard Model or a close extension. The original 
context of this is the weakly coupled heterotic string on a Calabi-Yau, where the 
matter spectrum is determined by the geometry of the manifold and the gauge 
bundles thereon. Realistic matter content requires substantial mathematical ef- 
fort; for an account of some recent developments [2] can be consulted. 

The style of top-down constructions changed with the discovery of D-branes 
and the realisation that Standard-like models could also be constructed within 
'intersecting brane worlds'. In these scenarios the Standard Model is located on 
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brane stacks that wrap cycles in the internal space, with the matter content de- 
termined by the intersection numbers of these cycles. Thus it is not just heterotic, 
but also (orientifolded) type II compactifications that can give pseudo-realistic 
physics. Such constructions have been developed both for toroidal orbifolds and 
smooth Calabi-Yaus, with a recent review being [3]. 

A general problem in top-down constructions is that of moduli stabilisation. 
The simplest string compactifications have very many moduli - massless scalar 
fields coupling with gravitational strength and determining the matter coupling 
constants. Such massless particles are however inconsistent with experiment and 
it is necessary to generate a potential for them. The difficulty is that moduli are 
associated with the geometry of the compactification, and the effects - preemi- 
nently fluxes - that stabilise the moduli also back-react on the geometry, making 
the resulting space difficult to study. 

Together with the brane world picture, this has led to 'bottom-up' approaches 
to string phenomenology [I]. As branes are localised, the field theory on them 
depends only on local geometry. A bottom-up physicist first builds a local Stan- 
dard Model, and only later worries about the global embedding. The limiting 
case of this is 'branes at singularities', where the entire low energy spectrum is 
determined by the nature of a pointlike singularity. In this context there have 
been recent attempts to find a Standard Model singularity [5]. A characteristic 
of this approach is that particle physics becomes an open string theory decoupled 
from the closed string dynamics associated with the compact geometry. 

Regardless of the philosophy, any realistic string compactification must even- 
tually account for the structure of the Standard Model. In this respect moduli 
stabilisation connects top-down and bottom-up constructions, as the moduli vevs 
determine the field theory coupling constants. One Standard Model coupling con- 
stant in particular need of explanation is the QCD 9 angle. This is of course a 
well-known problem with a well-known answer: a Peccei-Quinn axion [6]. There 
do exist other possibilities, reviewed in [TJ, but for this paper I shall simply assume 
the Peccei-Quinn solution to be correct. 

In the context of string theory, this creates a modulus anti-stabilisation prob- 
lem. There are many string theory axions that may in principle solve the strong 
CP problem. To do so, an axion must remain massless throughout the thicket 
of moduli stabilisation effects and down to the QCD scale. This problem has 
previously been brought up in [HI EJ- It is clean and sharply posed, as effects 
very weak on the Planck scale may be very large on the QCD scale. Given the 
necessity of moduli stabilisation, this question is best analysed within the context 
of string constructions for which all moduli have been stabilised. The purpose of 
this paper is to determine for such constructions the conditions under which a 
QCD axion will exist. 

The paper is organised as follows. I first review in section [2] the strong CP 
problem and the ways axions can arise in different string theory constructions. 
Section [3] investigates how to stabilise moduli while keeping axions sufficiently 
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light to solve the strong CP problem. I review some approaches to stabilising all 
moduli and show why their simplest versions contain no light axions. I investi- 
gate potentials with massless axions and show that in such cases there exist no 
super symmetric minima of the potential, but nonsupersymmetric minima may 
exist. I give conditions on the geometry for a QCD axion to exist, both at lead- 
ing and higher orders in the instanton expansion. Section H] addresses the axion 
decay constant. In the context of the exponentially large volume compactifica- 
tions of [Till [TT] I outline how a phenomenologically acceptable value of f a may 
be achieved together with the relationship f a ~ y/M S usYM P . 



2 Axions 

2.1 Axions and the Strong CP Problem 

In principle, the strong interactions can generate CP violation through an FF 
coupling, 



167T 2 _ 

However, observationally 9 is extremely small: \9\ < 1CT 8 . Given that at the 
above level, 9 is a coupling constant which can a priori take values anywhere 
between and 2tt, this seems unnatural. The strong CP problem is to explain 
this observation. 

There exist several proposed resolutions. In the context of string theory and 
string compactifications, the most natural is that due to Peccei and Quinn [6]. 
In this approach 9 is promoted to a dynamical field, with Lagrangian 

C = C SM + + Y^F^. (2) 

In (J2]) f a has dimensions of mass and is known as the axion decay constant. 
Conventionally a scalar has mass dimension one, and so we redefine a = 9f a . 
This gives 

C = C SM + \d,ad»a + —^—F^F^. (3) 

Z 107T J a 

In equation (J3]) there exists an anomalous global U(l) symmetry, a —>■ a + e. 
This symmetry is violated by QCD instanton effects, which break it to a discrete 
subgroup. These generate a potential for a, 

instanton ~ A QCD (l ~ cos (jj J ■ ( 4 ) 

In the absence of other effects, this potential is minimised at a = 0, setting the 
#-angle to zero dynamically. The mass scale for the axion a is 

A 2 

QCD / K x 

m a — . (5) 

Ja 
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(A more precise estimate replaces Aq CD by f^m^-, where f n ~ 90MeV is the pion 
decay constant.) The Peccei-Quinn symmetry is by its nature anomalous. For 
the axion to solve the strong CP problem, the leading anomalous contribution to 
the potential must be that of QCD instantons, otherwise the minimum will be at 
9^0. 

Phenomenologically, f a has only a narrow window of allowed values. The 
smaller the value of f a , the more strongly the axion couples to matter. The con- 
dition that supernovae cool by predominantly emitting energy through neutrinos 
(rather than axions) is equivalent to the constraint f a > 10 9 GeV. There is also a 
cosmological upper bound on f a . The axion field presumably starts its cosmo log- 
ical evolution with 9 at some arbitrary value between and 2n. Once the Hubble 
scale is comparable to the axion mass, the axion field oscillates and the energy 
density stored in the axion field is diluted with the expansion of the universe. The 
energy stored today thus depends on the axion mass, which is determined by f a . 
The requirement that axions do not overdose the universe leads to the constraint 
f a < 10 12 GeV, with such axions being potential dark matter candidates. 

The requirement f a < 10 12 GeV follows from the assumption of a standard 
cosmology. String compactifications generally have many moduli, which in super- 
symmetric scenarios typically have masses m ~ ITeV. These can be long-lived, 
causing problems with nucleosynthesis and giving low reheat temperatures. In su- 
persymmetric models, the axion always has a scalar partner (the saxion) with the 
associated cosmological problems. It has been argued (in particular see p^2l [13] ) 
that the cosmological problems with the scalars are more severe than those asso- 
ciated with the axion, and so the upper bound on f a should not be taken seriously 
without an associated resolution of the cosmological moduli problems. [T3] re- 
ports that models can be found in which both axion and saxion cosmological 
problems may be evaded with f a ~ 10 15 GeV. 

The cosmological history of the universe before nucleosynthesis is not known. 
However we are going to assume a standard cosmology and take the upper bound 
on f a seriously. This is mainly because low reheat temperatures (T ~ lOMeV) 
have a generic problem that at reheating the decaying modulus has an 0(1) 
branching fraction to gauginos and thus overproduces LSPs [TI[[T5]. 111 addition, 
low reheat temperatures make it hard to produce the observed baryon asymmetry. 
This suggests the universe should be hot at the weak scale, allowing a standard 
WIMP annihilation calculation and the possibility of electroweak baryogenesis, 
but also imposing the standard bounds on f a . Furthermore, this is not in obvious 
contradiction with the properties of the saxion - its matter coupling is also set 
by fa (rather than Mp), and low values of 10 9 GeV < f a < 10 12 GeV can easily 
give large reheating temperatures T ~ 10 5 GeV. 
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2.2 Axions in String Theory 



String compactifications generically contain fields a« which have ciiFF couplings 
and possess the anomalous global symmetry a —>■ a + e featuring in the axionic 
solution to the strong CP problem. We first enumerate possible axions, before 
considering their relation to the physics of moduli stabilisation. 



Axions in the Heterotic String 

In heterotic compactifications, the axions are traditionally divided into the uni- 
versal, or model- independent, axion and the model-dependent axions. The model- 
independent axion is the imaginary part of the dilaton multiplet, S = e~ 2 ^V ' + ia. 
It is the dual of the 2-form potential B 2 ^ v arising from the NS-NS 2-form field: 
da = e~ 2 ^ * dB^ v . The dilaton superfield is the tree-level gauge kinetic function 
for all gauge groups: 

C ~ Re{S) J F^ u F a ^ + lm{S) J F a ^F a ^ . 

Consequently in a realistic compactification there must always exist an clFqcdFqcd 
coupling. 

There are also the model independent axions, 6j, given by the imaginary parts 
of the Kahler moduli Tj. For a basis Sj of 2-cycles of the Calabi-Yau, = ti + ibi, 
with ti = f^. J the string frame volume of the cycle Sj and bi = L B 2 . At tree 
level, these have no couplings to QCD. However, such a coupling may be generated 
through the one loop correction to the gauge kinetic function. For the E s x E 8 
heterotic string, this correction is 

h = 3 + /3m, (6) 
h = S-m, (7) 

where 1 and 2 refer to the first and second E s respectively. The /3j are determined 
by the gauge bundles on the compactification manifold X. For gauge bundles V\ 
and V2, 

A = ^| e« A , (8) 

where is the 2-form associated with the cycle Sj. This can be derived by dimen- 
sional reduction of the Green-Schwarz term J B 2 A X 8 (Fi, F 2 ,1Z)\^ The axions 
associated to the Kahler moduli are called model-dependent as their appearance 
in / depends on the one-loop correction, which in turn depends on the specific 
properties of the compactification. 



2 Strictly, this only gives the correction to Im(/). The corresponding correction to Re(/) is 
however implied by holomorphy. 
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Axions in Intersecting Brane Worlds 



The discovery of D-branes [16] led to the extension of string model building 
beyond the heterotic string. The type II string theories, or more properly orien- 
tifolds thereof, can give rise to 'intersecting brane worlds'. In these the standard 
model is localised on a stack of branes while gravity propagates in the bulk. 
Light bifundamental matter arises from strings located at intersection loci and 
stretching between brane stacks. 

The bosonic action of a single Dp-brane is 

S P = - ~X +1 ( [ <F +1 Ze-+y/det{g + B + 2na>F) + i [ e B+2 ™' F A £ C q ) . 

(9) 

£ is the cycle wrapped by the brane and the sum is a formal sum over all RR 
potentials in which only relevant terms are picked out. 

The kinetic term F fll/ F IJ,u comes from the DBI action and the instanton action 
F A F from the Chern-Simons term. The gauge coupling corresponds to the 
inverse volume of £ and the 6 angle to the component of C p _3 along S. These 
fields pair up to become the scalar component of the chiral multiplet which is the 
gauge kinetic function of the resulting gauge theory. 

As IIB compactifications are our main focus, we shall be more explicit here. In 
principle, IIB string theory allows, consistent with supersymmetry, space-filling 
D3, D5, D7 and D9-branes. However, in an orientifold setting we are restricted to 
either D3/D7 or D5/D9 pairs. We shall interest ourselves in the former case. The 
branes must wrap appropriate cycles to cancel the negative charge and tension 
carried by the orientifold planes; we assume this has been done. 

In D3/D7 compactifications, the relevant superfields are those of the dilaton 
and Kahler moduli multiplets. Their scalar components ar^f] 

S = e^ + ico, (10) 
Ti = n + ici. (11) 

Co is the Ramond-Ramond 0-form and = g s the string coupling. For E$ a 
4- cycle of the Calabi-Yau, 




where l s = 2ny/a' denotes the string length. Note that the Kahler modulus 
Ti involves the Einstein, rather than string, frame volume of a 4-cycle. This is 

3 Technically, this is only for the case that ft, I' 1 = 0. This will not be important for the issues 
we discuss, and so we will use this simplifying assumption. The correct expressions under more 
general circumstances can be found in p~7l [18] . 
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most simply understood as the requirement that the gauge kinetic function be 
holomorphic in the chiral superfields. Indeed, S is the universal gauge kinetic 
function for D3-branes, whereas Tj is the gauge kinetic function for the field 
theory on a D7-brane stack wrapping the 4-cycle Sj. 

The axionic couplings arise from the Chern-Simons term in the action. For 
D3-branes, this gives 



S 



while for D7-branes 



FF 



FF 



Ci_ 

An 



F A F. 



(12) 



(13) 



By expanding the DBI action, we obtain the field theory couplings 



D3 



27 



and — 



D7 



2tt' 



There exists a similar story for IIA intersecting brane worlds, where the Stan- 
dard Model must be realised on wrapped D6-branes (a Calabi-Yau has no 1- or 
5-cycles to wrap D4- or D8-branes on). The gauge coupling now comes from the 
calibration form Re(f2) and the axion from the reduction of the 3-form potential 

C 3 , 

= / Re(ft) and 6 
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Our interest is in the interaction of axions with moduli stabilisation and super- 
symmetry breaking, to which we now turn. 



3 Axions and Moduli Stabilisation 

It is obvious from the above that potential axions are easily found in string com- 
pactifications; indeed, they are superabundant. For axions to solve the strong 
CP problem, they must also be light, with QCD instantons giving the dominant 
contribution to their potential. Light axions are not in themselves problem- 
atic. Pure type II Calabi-Yau compactifications have many axions, which remain 
exactly massless as a consequence of four-dimensional Af = 2 supersymmetry. 
However, the same Af = 2 supersymmetry that guarantees the axions remain 
massless also guarantees a non-chiral spectrum with the axions' scalar partners 
massless. These are modes of the graviton and will lead to unobserved fifth forces. 

More realistic string constructions have Af = 1 supersymmetry in four di- 
mensions, allowing a potential to be generated for the moduli. To avoid bounds 
from fifth force experiments, the size moduli - the saxions - must at a minimum 
receive masses at a scale ttit > (100/rni) -1 . However, the expected scale is much 
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larger: typical constructions give moduli masses comparable to the supersymme- 
try breaking scale, and the cosmological moduli problem [19j 120] suggests that 
in fact rriT > lOTeV. In recent years there has been much progress in moduli 
stabilisation [21], [22], with physics such as fluxes and instantons used to lift the 
degeneracies associated with the geometric moduli. In the context of strong CP, 
this same progress creates a modulus anti-stabilisation problem. There are many 
stringy effects that can generate a potential for a putative QCD axion. These in- 
clude worldsheet instantons, D-instantons and gaugino condensation, all of which 
are often invoked to stabilise the geometric moduli present. If any one of these 
effects is more important for a given axion than QCD instantons, that axion does 
not solve the strong CP problem. 

Axion potentials generally come from nonperturbative effects whose magni- 
tude is exponentially sensitive to the values of the stabilised moduli. The analysis 
of such effects is therefore best performed in a context within which all moduli 
are stabilised, and we first review mechanisms to achieve this. The most devel- 
oped scenarios are those within IIB flux compactifications. We will work mostly 
with these, although we will along the way obtain a no-go theorem applicable 
to all string compactifications. IIB compactifications have the advantage that 
the stabilisation of dilaton and complex structure moduli can be studied ten- 
dimensionally [21], while the back- reaction of the fluxes is relatively mild: the 
internal geometry simply becomes conformally Calabi-Yau. 



3.1 Review of Moduli Stabilisation 

We shall orient our subsequent discussion around two approaches to stabilising 
all moduli, the well-known KKLT scenario [22] and the exponentially large vol- 
ume models of [TO], [IT]. The point of departure is the flux compactifications of 
Giddings, Kachru and Polchinski [21]. These are compactifications of a IIB ori- 
entifolqj in the presence of localised sources (D3/D7 branes and 03/07 planes) 
with non- vanishing 3- form flux G3 = F3 — SH3. Here F3 and if 3 are RR and NSNS 
3-form fluxes with S the dilaton-axion. The fluxes generate a superpotential[23] 

W = ^Jg 3 AQ, (14) 

with Q the unique holomorphic (3,0) form and the Kahler potential given by 

A -21n(V)-ln( [qaQ ) -\n(S + S). (15) 



M 2 P 



This Kahler potential is no-scale for the Kahler moduli. The dilaton and complex 
structure moduli appear in the superpotential and are stabilised at leading order 



or more generally F-theory. 
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in the g s and a 1 expansions. Their masses scale as 

N 

m s ~ ~ ~ftE Mp ' ( 16 ) 

where N is a measure of the number of units of 3- form flux and R is the Calabi- 
Yau radius in units of l s . 

At this level, the no-scale property implies the Kahler moduli remain unsta- 
bilised. The dilaton and complex structure moduli are integrated out to focus 
on an effective theory for the Kahler moduli. Although absent at tree level, the 
Kahler moduli can appear nonperturbatively in the superpotential through brane 
instantons [2l] or gaugino condensation, 

w = w Q + J2 A i e ~ aiTi > ( 17 ) 

i 

where Oj = 2tc(^-) for brane instantons (gaugino condensation). The KKLT 
proposal is to stabilise the Kahler moduli by solving D^W = for all i. The 
resulting 4- cycle sizes are 

1 , fw 



^ ln UJ- (18) 

The cycle size goes logarathimically with Wq, and so to trust the supergrav- 
ity approximation Wq must be extremely small. We may usefully rephrase this 
condition. The KKLT construction relies on nonperturbative corrections to the 
scalar potential dominating the perturbative corrections. This is equivalent to 
the requirement that the 'correction' to the superpotential is comparable to, or 
larger, than the tree-level term. In order for this to hold at (moderately) large 
T, Wo must be very small, which is to be achieved by tuning fluxes. 

However, for almost all values of Wo ('almost' can be made precise as in [25] ) 
there exists no reliable regime in which this requirement is met. Perturbative 
corrections, originating from corrections to the Kahler potential, are then essen- 
tial for the study of the moduli potential. It was shown in [TUJ [11] that the 
incorporation of an a' 3 correction to the Kahler potential [26] , 

£ = -21n(V)->/C = -21n(V + -M , (19) 

V 2glJ 

combined with the same nonperturbative superpotential corrections 

W = W + ^2A i e- ttiTi i 

i=2 

leads, subject to a necessary condition h 2,1 > h 1,1 > 1, to a minimum of the 
potential at exponentially large volumes: 

c 

V ~ Woe 3« , 
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for a model-dependent constant c. This and other properties of the minimum 
follow from an explicit study of the scalar potential, but a summary of the results 
is as follows: 

1. The moduli divide into one large modulus Tj, and h 1 ' 1 — 1 small 'blowup' 

2 

moduli Tj. Whereas the former is exponentially large, T\, ~ V^, the latter all 

3 

have Tj ~ ln(V). The origin of this is as follows. If we take r b 2 ~ V ^> 1 and 
fixed, the other moduli Tj minimise their potential at D^.W = + O (i). 
Neglecting numerical factors, this generates an effective potential for V of 

where the ^ comes from the Kahler correction. V is then dynamically 
stabilised at V 3> 1. 

2. The stabilised volume is exponentially sensitive to the stabilised string cou- 
pling, allowing a natural generation of hierarchies: the structure of the 
potential is such as to create a hierarchical separation of the string and 
Planck scales. 

The Kahler correction (|T9l) - of course only one of many - plays an essential role 
in this construction. The effects of other possible Kahler corrections, such as from 
warping effects or other IIB 0(a' 3 ) corrections, have been extensively discussed 
in [11] (see also [28J). We will comment further on these in the appendix. Here we 
simply note that the only hope of incorporating Kahler corrections in a controlled 
manner is via a small expansion parameter. In this case V>1 and the expansion 
parameter is y. Of known corrections, (|T9l) gives the leading contribution to the 
scalar potential in the h expansion. 

Finally let us note that the two constructions above both give rise to AdS 
minima, which are respectively supersymmetric and non-supersymetric. It is 
necessary to include a phenomenological uplift term, arising from e.g. a D3-brane 
in a warped throat, in order to match the observed cosmo logical constant. In the 
KKLT scenario, it is the uplift that gives rise to soft supersymmetry breaking 
terms. For the exponentially large volume compactifications, the supersymmetry 
breaking is inherited from the no-scale structure and the effects of the uplift are 
sub dominant [TT1 129] . 

3.2 The Simplest Scenarios: Why Axions are Heavy 

We now turn to the axions. Our first point is that the simplest versions of the 
above scenarios lack a QCD axion. All potential axions receive a high scale mass 
and thus cannot solve the strong CP problem. For simplicity we concentrate on 
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the KKLT construction, but a very similar argument holds for the exponentially 
large volume compactifications. 

We start by asking whether QCD is to be realised on D3 or D7 branes. If 
we were to use D3-branes, the QCD axion would be the imaginary component of 
the dilaton multiplet, S = + ico. However, as indicated in ({IB]) this multiplet 
is stabilised at tree-level by the fluxes, with a mass iris ~ 

^M P , with R the 

radius in units of l s . This tree-level stabilisation may seem at odds with the 
axionic shift symmetry cq — > Cq + 2n. However, the shift symmetry is a subgroup 
of the fundamental SL(2, Z) symmetry, under which the fluxes also transform. 
By nature a duality, this is invisible in the low-energy theory and cannot protect 
the axion from acquiring a mass. As R < 5 in KKLT, the axion obtains a mass 
m Co > 10 15 GeV and cannot be a QCD axion. For the exponentially large volume 
compactifications, R is larger but the conclusion unchanged: QCD on a D3 brane 
stack is inconsistent with the existence of a Peccei-Quinn axion. 

This implies that QCD ought to be realised on D7 branes. The axions are 
now the imaginary parts of the Kahler moduli, and the instanton effects used to 
stabilise these moduli will also give the axions a mass. To estimate the scale of 
this mass, it is simplest just to construct the potential explicitly. 

As above, we take the superpotential 

W = Wo + J2 A i e M-aiTi), (20) 

8=1 

with Kahler potential 

/C = -21ogV. (21) 

Note that JC = /C(Tj + T«), and so the axions do not appear either in JC or its 
derivatives. The supergravity F-term potential is 

V = e K (Kf j D i WD 3 W - 3|H/| 2 ). (22) 

The no-scale property of the Kahler potential simplifies (1221) to 

V = e K ^d.WdjW + KP ((dJQWdjW + (d-^WdjW)) . (23) 

It is a property of the Kahler potential (T2T1) that JC^diJC = —2tj. f[2"3"j) becomes 

V = e K (^Kf j a i a j {A i A j e- a ^- a ^ + A i A j e- aiTi - aiTj> ) 

-2aiTi (wAie-^ + WA^-^ J . (24) 

It is easy to extract the axionic dependence of the potential ( 1241) . JC and its deriva- 
tives are all real and phases only come from the superpotential. The potential 



12 



becomes 



V = e K (jC i] (2a i a j \A i A j \e- am - a ^ cos(a i 6' i + aj 0j + 7y)) (25) 
—AaiTi\W Ai\e~ aiTi cos(a i 6 l j + — Aa^AiA^ cos(cij0j + + 7^) J . 



9i denote the axions and the phases 7^ and $ come from the phases of AiAj and 
AiW respectively. The axionic mass matrix is 



M « = mk/ (26) 



and we obtain^! 

Mfj ~ Oiaj-y m< „, (27) 

with V^.j n the magnitude of the potential at the AdS minimum. In KKLT, 

D T W = for all % and V min 3^. As there are ft, 1 ' 1 independent phases, 

there is no reason for M?- to be degenerate and we expect all eigenvalues to be 
0(a 2 V m i n ), where a is the typical magnitude of the e^. 

The determination of physical masses also requires the Kahler potential. In 
general there is no explicit expression for the overall volume V in terms of 4-cycle 
volumes Tj. The Kahler metric may however be written as [32] 



G.- 1 3 fk -,v k ?,h f m t n h- t p t q 



If V ~ (a few)/f, /Cjj ~ 0(1) and the mass matrix M?- gives a good estimate of 
the scale of axion masses. If V ^ if, then /Qj <C 0(1), and M?- underestimates 
the axion masses. Mf- could only overestimate the axion masses if V 1. This 
realm of moduli space is not accessible in a controlled fashion and we do not 
concern ourselves with it. 

In units where Mp = 1, we therefore have 

m n ~ m Ci ~ diVKun ~ • (29) 

The axion masses are consequently set by the value of the tree-level superpotential 
Wo- This also determines the vacuum energy and, implicitly, the energy scale of 
supersymmetry breaking required to cancel the vacuum energy. TeV-scale soft 
terms require hierarchically small Wo- For the (gravity- mediated) case studied in 
[30] . this required W ~ 10~ 13 , with 

m Ti ~ m a% ~ m 3 / 2 ~ lOTeV. (30) 



5 We emphasise this does not mean the mass matrix is M 2 = a® a: simply that terms may 
receive an enhance by factors appearing in the exponential. 
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This scale is vastly greater than that associated with QCD instanton effects, and 
thus the axions are incapable of solving the strong CP problem. We could insist 
on a QCD axion, and require that Wo be sufficiently small that QCD instantons 
dominate over the D-instanton effects of moduli stabilisation. This would require 
Wq ~ 10 -40 . However, this scenario is entirely excluded as the size moduli are 
light enough to violate fifth force experiments and the susy breaking scale would 
be (9(10~ 14 eV). Consequently, in the simplest KKLT scenario it is impossible 
to generate a QCD axion. The D3 axion receives a high scale mass from fluxes, 
whereas the instanton effects give the D7 axions large masses comparable to the 
size moduli. 

A similar argument holds for the exponentially large volume compactifica- 
tions. As in KKLT, the D3 axion receives a large flux-induced mass. The 'small' 
cycles are stabilised by instanton effects, and these give the corresponding axions 
masses of a similar scale to the size moduli, m ai ~ m Ti ~ m 3 / 2 . One difference is 
that there is a modulus, the 'large' modulus r b , which need not appear in the su- 
perpotential. It is stabilised through the Kahler potential, and while it is massive 
its axionic partner indeed remains massless. However, this cycle is exponentially 
large, and any gauge group supported on this cycle is far too weakly coupled to 
be QCD. The same conclusion holds: the simplest version of this scenario does 
not generate a QCD axion. 

The above formulates the 'modulus anti-stabilisation problem': naive scenar- 
ios of moduli stabilisation are incompatible with a QCD axion. 

We next examine an apparent solution to this problem that in fact has a subtle 
flaw. We want a way to stabilise moduli without stabilising the axions. Axions 
correspond to phases in the superpotential and do not appear in the Kahler 
potential. If we included a multi-exponential term e~ atTi in the superpotential, a 
massless axion would certainly survive, as at least one phase would be absent. As 
the size moduli all appear in the Kahler potential, by solving the F-term equations 
we may hope to stabilise the size moduli while leaving the axions massless. 

To illustrate this idea, let us consider a toy model, 



The Kahler potential (I3TT) is that appropriate for toroidal orbifolds, with V = 
ti^2^3- To understand where the superpotential could arise from, we can hypoth- 
esise that the cycle (1+2+3) is the smallest cycle with only two fermionic zero 
modes, and that instantons wrapping (for example) the cycle (2+3) all have more 
than two zero modes and do not appear in the superpotential. However, we are 
not here really concerned with the microscopic origin of the superpotential: at 
this level we simply regard equations (J3~lj) and fl32l as defining the model. 
The F-term equations D Tl W = D T2 W = D Ta W = give 



W 



- ln(Tx + 7\) - ln(T 2 + T 2 ) - ln(T 3 + f 3 ) 

W + Ae- 2 ^ +T2+T *\ 



(31) 
(32) 



- 27r^e- 27r(Tl+T2+T3) - 



1 



(Wo + Ae~ 2 ^ +T ^) 



(33) 







T 1 + T 1 
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-±-^(W + Ae- 2 ^ +T ^) = 0, (34) 

2 + ^2 

-l^(W + Ae- 2 ^ +T2+T ^) = 0. (35) 

3 + T 3 

These immediately imply t\ = r 2 = r 3 and equations (1331 to fl35|) collapse to 

2nAe~ 67TT1 e- 2m{ei+e2+6:i) + — (W + Ae~^ Tl e- 2m{(>1+e2+ez) ) = 0. (36) 

While the sum 0i + 02 + 03 is fixed, there are clearly two axionic directions not 
relevant for the solution of the F-term equations. On the other hand, there is 
a unique value for the size moduli such that the F-term equations are solved. 
Except for the massless axionic directions, the scales of the masses are unaltered 
from above, and we would expect 

W ° n (i>7\ 

m Ti ~ m ei+ g 2+ e a ~ — , m ei „ 02 = m 8l ^e :j = 0. (37) 

As this is supergravity rather than rigid supersymmetry, there is no contradic- 
tion in having a mass splitting for the multiplet in the presence of unbroken 
sup er sy mmet ry. 

While this seems promising, there is in fact a serious problem with the 
above. Even though all F-term equations can be solved, numerical investiga- 
tion shows that at the supersymmetric locus the resulting scalar potential is 
tachyonic, with signature (+, — , — ). Although supersymmetry ensures the mod- 
uli are Breitenlohner-Freedman stable [31J, this notion of AdS stability ceases to 
be relevant after the (necessary) uplift. 

We now show that these tachyons are in fact generic for any attempt to 
stabilise the moduli supersymmetrically while preserving unfixed axions. 



-2nAe- 27T( - Tl+T2+T3 ~> - 
-2nAe- 2 ^ Tl+T2+T3) - 



3.3 A No- Go Theorem 

We suppose we have an arbitrary M = 1 supergravity with moduli fields, $ a , 
Tp = Tp + ibp, where the bp are the axions. We write the superpotential and 
Kahler potential as 

W = W(<5> a ,Tp), (38) 
K = !C($ a ,Tp + Tp). (39) 

The Peccei-Quinn symmetry bp — > bp + ep implies the form of fl3"9~l) should hold 
in perturbation theory. 

We further suppose we have solved 

D$ a W = and D T0 W = (40) 
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for all a and /?, but that at least one axion b u = is unfixed: the solution 

to (I4TJ]) is independent of (b u ). 

We redefine the basis of chiral superfields so that there exists a superfield T u 
with b u = Im(T M ), 

T\ — > T u , 
T 2 - T 2 , 

T„ - T„. (41) 

This is a good redefinition as it does not affect holomorphy properties. 

As the solution to all F-term equations is independent of b u , b u is a flat 
direction of the potential (1221) at the supersymmetric locust The potential at 
the supersymmetric locus is given by 

V = -3e lc \W\ 2 . (42) 

As b u does not appear in /C, it follows that if b u is a flat direction \W\ must be 
independent of b u . Up to one exception this then implies that W is independent 
of b u . 

The sole exception is if b u purely represents an overall phase, i.e. W = 
e -aT u w j^] 1 no constant term. In a IIB context, this may potentially arise if 
the flux superpotential exactly vanishes due to a discrete symmetry (32], while 
a combination of non-perturbative effects and Kahler corrections stabilise the 
Kahler moduli. While potentially interesting, this is an exceptional case and we 
do not analyse it further. 

If W has no explicit dependence on b u = Im(T u ), it follows by holomorphy 
that it also has no explicit dependence on t u = ~Re(T u ) and hence on T u . Therefore 

d Tu W = 0. (43) 

However, as Dt u W = 0, it follows that at the supersymmetric locus, 

either {d Tu K) =0 or W = 0. 

The latter is overdetermined and non-generic, so we first focus on (8t u K) = 0. 

Direct calculation now shows that the t u direction is tachyonic at the super- 
symmetric locus. To see this, note that from the scalar potential (|2"2]1 

d Tu V = e K K, fj (d Tu {DiW)D 3 W + DiWd Tu (DjW)) - 3(d Tu K)e K WW. (44) 

6 The requirement of flatness is stronger than the requirement that the axion simply be 
massless. Flatness is the right requirement, as if an axion is fixed in any way it does not solve 
strong CP. 
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We have used d Tu W = and have only kept terms that will give non- vanishing 
contributions to d Tu d Tu V at the supersymmetric locus. Expanding DiW and 
again using d Tu W = 0, ( 1441) simplifies to 

d Tu V = e K Kfi (dr v (diK)W{DjW) + D.Wd^iqw) - 3(d Tu )C)e K WW. (45) 

If we again only keep terms non- vanishing at the supersymmetric locus, the second 
derivative is 

d Tu d Tu V = e K Kfi (2d Tu {d i K)d Tu {d- j K)WW) - 3(d Tu d Tu JC)e K WW. (46) 
Now, as t u = \{T U + f u ), 

d Tu JC(T + T) = 2d Tu JC(T + T), 

and we have 

d Tu d Tu V = 4e !C WW(2IC^IC tu IC u] -3IC U u) 

= -Ae^WWIC^, (47) 

where we have used = Kiy As is a metric, K U u is positive definite and it 
follows that the r u direction is tachyonic. 

Now consider the W = case. As indicated above, this is non-generic: even 
if W originally vanishes, it is expected to receive non-perturbative corrections 
which make it non-vanishing. Even so, it follows easily that 

d Tu d Tu V = 0, (48) 

and so the t u size modulus is massless, leading to unobserved fifth forces. 

The above gives a no-go theorem: there does not exist any supersymmetric 
minimum of the F-term potential consistent with stabilised moduli and unfixed 
axions. 

It is in the nature of no-go theorems that they admit loopholes, so let us 
discuss ways around this result. One point to consider is the form of the Kahler 
potential, as we have used in the above argument the fact that 

/C = /C(T + f). (49) 

While true in perturbation theory because of the axionic shift symmetry, this 
equation will break down nonperturbatively, and the argument showing that the 
r u direction is tachyonic will no longer hold. However, this same breakdown will 
cause /C, and hence the potential V, to depend on the axion. As this lifts the 
required axionic flat direction, the no-go theorem will cease to apply. 

A second loophole is that although the F-term potential might be tachyonic, 
the D-term potential might come to the rescue. For example, a Fayet-Iliopoulos 
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term might have exactly the right structure to render the supersymmetric locus 
an actual minimum of the full potential. However this seems implausible in the 
presence of many tachyonic directions. A similar approach would be to try and 
set W = and then rely entirely on D-terms to stabilise the moduli. Another 
possibility (discussed recently in [33]) is that an anomalous U(l) might remove the 
tachyonic directions. While the massive gauge boson will eat the axionic degree 
of freedom, an axionic direction may survive in the phase of a scalar charged 
under the U(l). 

A third loophole is that stability might not rely on the existence of an actual 
minimum for the potential. Equation ( H7|) involves the Kahler metric /C u „. The 
kinetic term for t u is KL U ud^T u d^r u . If we just consider the t u direction, the 
physical mass is therefore 



where vhbf is the relevant Breitenlohner-Freedman bound. As tachyonic modes 
can be stable in AdS, one could argue that it is sufficient simply to solve the 
F-term equations and not to worry about whether the resulting locus is an actual 
minimum of the potential. 

While this point is more substantial, it does not resolve the problem. The real 
world is not AdS, and for stability requires a positive definite mass matrix. For 
any realistic model, the vacuum energy must be uplifted such that it vanishes. 
After this uplift, the extra geometric advantages of AdS go away and the tachyons 
can no longer be supported. As there may be many tachyons present - one for 
each massless axion - the entire problem of moduli stabilisation must necessarily 
be solved over again in the uplifting. As the uplift is generally the least controlled 
part of the procedure, this seems a hard problem. While the uplift may remove 
the tachyons, at the present level of understanding this seems pure hypothesis. 
It is then very unclear how useful the original supersymmetric AdS saddle point 
is, and whether it is a suitable locus to uplift. 

The argument above suggests that supersymmetric solutions are unpromising 
starting points from which to address the strong CP problem. Either all moduli 
appear in the superpotential, in which case there is no light axion, or a modulus 
is absent from the superpotential, in which case the potential is tachyonic. The 
fourth and most obvious loophole is then to give up on the requirement of super- 
symmetric minima, and search for nonsupersymmetric minima of the potential 
with massless axions. 

We shall consider this point in the next section. Before we do so we also 
observe that while our focus here is on axions, the argument above is also relevant 
for moduli stabilisation. For example, in the weakly coupled heterotic string, 
the one-loop corrections to the gauge kinetic function dSJ) imply that gaugino 
condensation generates a superpotential 




-2e K WW 




(50) 



W r 



n.p. 



(51) 
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It has been proposed to use the superpotential (I5"T|) . together with a constant term 
Wq, to stabilise the dilaton and Kahler moduli by solving D S W = D Tl W = 0. 
However, there is only one phase - and hence only one axion - explicitly present in 
( I5T1) . The above argument shows that the resulting scalar potential will actually 
be tachyonic at the supersymmetric locus, with signature (+, —,...,—). 

A similar result will apply to the recent study of IIA flux compactifications 
with all NSNS and RR fluxes turned on. In this context it is also found that the 
solution of the F-term equations is independent of many of the axions present 
(those associated with the C3 field). The above implies that as long as h 2,1 ^ 
the supersymmetric locus is always tachyonic, with one tachyon present for every 
massless axion. Tachyons have been recognised in particular models [33, 34J , but 
from the above they would seem to be very generic. 

3.4 Non-supersymmetric Minima with Massless Axions 

The above no-go theorem shows that supersymmetric moduli stabilisation is not 
a good starting point from which to solve the strong CP problem. This implies we 
ought to consider non-supersymmetric moduli stabilisation. That there is no no- 
go theorem for non-supersymmetric minima with massless axions can be shown 
by construction: for example, the exponentially large volume compactifications of 
[TOf [TT] all contain a massless axion associated with the large cycle controlling the 
overall volume. As indicated above, this cannot be a QCD axion, as any brane on 
this cycle is very weakly coupled. If we want to try and force this cycle into being 
a QCD axion, we can tune the parameters to force the minimum of this potential 
to relatively small volumes. That this is possible can be confirmed numerically. 
Another possibility would be a purely perturbative stabilisation of the volume 
modulus, solely using Kahler corrections (in which axions do not appear). This 
has been discussed in [23 EE], although without an explicit example. 

We shall not dwell on these possiblities. First, because the resulting axion 
decay constant would be, as we shall see in section [U close to the Planck scale 
and outside the allowed window and secondly, because at such small volumes 
there is no good control parameter. As we require non-supersymmetric minima 
we shall base our discussion around the large- volume models of [TOj, [H] . For now 
we only discuss moduli stabilisation but in section H] we shall show that these can 
also realise phenomenological values for /„. 

We illustrate the discussion with a three-modulus toy model, in which we 
assume the volume may be expressed in terms of 4-cycles as 

V = (T x + Ti)i - (T 2 + T 2 )l - (T 3 + t 3 )l. (52) 

We use three moduli as this turns out to be the minimal number required for our 
purposes: clearly, this is no significant restriction. Expressed in terms of 2-cycles, 
( 152"j) corresponds to 

V = X(tl-t 3 2 -4). (53) 
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We note (153]) satisfies the requirement that ^ > ^ t . have signature (+,—,—). We 

may perhaps think of this toy model as a P 3 with two points blown up. Denoting 
the cycles by 1, 2 and 3, the Kahler potential is 

K = -2 In ( (T a + f x )§ - (T 2 + T 2 )§ - (T 3 + T 3 )§) - -g|-, (54) 

where we have included the leading large-volume behaviour of the a' 3 correction 
of g s is fixed by the fluxes and in (1531) should be regarded as a tunable 

parameter. For superpotential, we shall take 

W = W + e-^ {T2+n) . (55) 

This could arise from gaugino condensation on a stack of n branes wrapping the 
combined cycle 2+30 QCD will be realised as a stack of branes wrapping cycle 
3. 

In the limit V ^> 1 with r 2 and r 3 small, the leading functional form of the 
scalar potential is (omitting numerical factors) 

T/ (y/^+y/^)e-^ +T3) (r 2 + r 3 )e-^+^ £ 

^ = ^1 1 — ^ — • (56) 

The minus sign in ( 1561) arises from minimising the potential for the axion Im(T 2 + 
T 3 ). The axions Im(Ti) and Im(T 2 — T 3 ) do not appear in (156]) and are unfixed. 
By considering the limit V — > 00, 27r(T2+T3) ~ In V, it follows that as V — > 00 the 
potential (156]) goes to zero from below. As by adjusting g s we can make the third 
term of (I56p arbitrarily large, we can ensure the potential remains positive until 
arbitrarily large volumes, and thus any minimum will be at exponentially large 
volumes. 

Is there a minimum? The potential is clearly symmetric under r 2 <-> r 3 , and 
the potential restricted to the locus r 2 = r 3 indeed has a minimum at exponen- 
tially large volumes. Because of the symmetry r 2 <-> r 3 , this 'minimum' is also 
a critical point of the full potential. However, it is not a minimum of the full 
potential. At fixed r 2 + r 3 and fixed V, (15B1) depends only on ■ s /tq, + For fixed 
r 2 + T 3, this is maximised at r 2 = r 3 , and the mode r 2 — r 3 is tachyonic at this 
locus. We have not investigated whether this tachyon satisfies the Breitenlohner- 
Freedman bound for AdS stability. This is for the same reasons as above: once 
we uplift, the geometric protection of AdS ceases to be relevant. Consequently 
the fields in (1561) run away either to r 2 = or r 3 = 0, where one of the blow-up 
cycles collapses. 

This result shows that the above toy model does not, by itself, have a minimum 
of the potential with a massless QCD axion. We may ask whether this is a feature 



7 The use of gaugino condensation rather than instanton effects is necessary to ensure that 
the cycle 2+3 is large enough to contain QCD. 
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of the geometric details of the model - for example, whether a different choice 
of triple intersection form in (|53|) would alter this result. We have investigated 
several other toy models without finding a minimum, and while we have no proof 
we suspect none exists so long as the Kahler potential is given by ( 1541) . 

This is bad news, but it is controllable bad news. The instability above is 
very particular: there is no instability either for the overall volume or for the sum 
of the blow-up volumes T\ + r 2 , but only for the difference T\ — r 2 . The effect 
of the instability is to drive one of the blow-up cycles to collapse. Consequently, 
the instability can be cured by any effect that becomes important at small cycle 
volume and prevents collapse. 

For example, the addition of a term 

1 + ^ 



v^V 3 ^V 3 

to fl56l) would obviously stabilise the cycles r 2 and r 3 against collapse and generate 
a minimum of the potential. As this term does not affect the argument that in the 
V —>■ oo limit the potential approaches zero from below, the resulting minimum 
would still be at exponentially large volume. 

We discuss in greater detail in the appendix which Kahler corrections are and 
are not allowed. For now we note that terms of the form (1571) may be generated 
from a correction to the Kahler potential, 



K, + 5K. = -2 ln(V) + ^ + (58) 

For simplicity we have kept the r 2 <-> T3 symmetry. Such a correction is motivated 
by the fact that it gives corrections to the Kahler metrics /C 22 and /C33 suppressed 
by factors of g 2 for the field theory on the relevant cycle. More specifically, 



/C 2 o + <5/C 2 2 — — ; — — I 1 = — I . 

2^2F 2 V V 12 



(59) 



As r 2 = p for a brane wrapping the cycle 2, the correction is suppressed by g 2 . 

The inverse metric involves an infinite series of terms diverging in the r 2 — > 
and T3 — > limit. For example, 



, = 2jW( e ^_ y (60) 

3 V 12V2t 2 288rf 1 

This is easy to understand: at r 2 = ^^75? ^ ne Kahler metric /C 22 goes to zero and 
the inverse metric diverges. This divergence can be seen by resumming (|60l) . In 
the physical potential, this divergence will create a positive wall at finite values 
of r 2 and r 3 . The positivity can be seen from the fact that the divergence in /C _1 
will only appear in the term 

e^VWiWDjW, (61) 
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which is manifestly positive definite. Consequently the potential will diverge 
positively at finite values of r 2 and r 3 , and so a stable minimum must exist for 
both r 2 and r 3 . 

We have now outlined, in the context of the scenario of pm EL a wa Y to 
stabilise all the size moduli while keeping a massless QCD axion left over. We 
would like cycle 3 to support QCD: by adjusting e, we can always make the 
correction fl58|) sufficiently large to ensure that t 3 is stabilised with the correct 
size for QCD. For intermediate string scales, this requires r 3 ~ 10. As the 
actual correction would be very hard to calculate, at this level we just adjust e 
phenomenologically. Of course, the complexity of a real model is much greater 
than that of (1581) . However we note again that, even though the corrections 
cannot be calculated, our proposal for moduli stabilisation only requires that 
they exist and come with the right sign to prevent collapse. 

While the above has been with a toy example, the above approach will apply 
to any model in which the moduli are stabilised along the lines of [T0~l[TT] . Keeping 
an axion massless introduces an instability causing a blow-up cycle to want to 
collapse. Kahler corrections that become important at small volume can stabilise 
this cycle but will not affect the overall structure of the potential, and in particular 
will not affect the stabilisation of the volume at V 3> 1. 



3.5 Higher Instanton Effects in the Axion Potential 

We have given above a Kahler potential and superpotential that will stabilise the 
moduli while containing a candidate QCD axion. The nonperturbative terms in 
the superpotential are in general just the leading terms in an instanton expansion. 
Even though the higher order terms may be highly suppressed and irrelevant to 
moduli stabilisation, they could still lift the flat direction associated with the 
massless axion. Given that Aq C d -C M p , even highly sub leading terms could 
dominate over QCD instantons. 

Let us estimate the general magnitude of such instanton effects. The mag- 
nitude of brane instantons depends on the volumes of the cycles they can wrap. 
Generally there will be many such cycles, whose sizes depend on the stabilised 
moduli, but minimally there must always exist the cycle which support the QCD 
stack. It follows from the DBI action that the gauge coupling for a D7-brane 
stack is^| 

1 Re(T) . 4tt „ /m , 

- = — U. = Re T = r. 

9 4tt g 2 

This defines the gauge coupling at the high scale where the effective field theory 
becomes valid. When m s ~ Mp, this is in essence the string scale, but if m s <C 



8 v3: This corrects a factor of 2 from earlier versions. To deduce the gauge coupling from 
reduction of the DBI action, it is necessary to ensure the gauge group generators are normalised 



according to the phenomenology conventions Tr{T a T 13 ) — ^S a 



22 



Table 1: Cycle sizes and instanton amplitudes for various UV scales 



Euv 


10 8 GeV 


10 10 GeV 


10 12 GeV 


10 14 GeV 


10 16 GeV 


a QCD^Vv) 


15.8 


18.1 


20.4 


22.7 


25 


Re(T) = a' 1 


15.8 


18.1 


20.4 


22.7 


25 




7.7 x 10" 44 


2 x 10" 50 


2 x 10" 56 


6 x 10" 63 


6 x 10" 69 




5.9 x 1(T 87 


4 x lO" 100 


4 x 10" 112 


3.6 x 10" 127 


3.6 x 10~ 139 



Mp, the difference between uikk and m s become significant. It is a subtle issue 
whether m s or m KK is the appropriate high scale. If QCD is supported on a small 
cycle within a large internal space, the KK modes associated with the bulk will 
be uncharged under QCD and will not contribute to the running coupling. KK 
modes of the QCD cycle will contribute, but these will be at masses comparable 
to the string scale. In considering the running coupling, we therefore use m s as 
the high scale rather than ttikk- We consider a wide range of string scales and 
take a sampling of high scale values from 10 8 — > 10 16 GeV. Given a string scale, 
the internal volume is determined by m s ~ 

The QCD coupling runs logarithmically with energy scale, with 

aQ^(10 2 GeV) ~ 9 and a^ D (10 16 GeV) ~ 25. 

The required high-scale couplings and cycle sizes are given in table (TJ together 
with the action for a D3-brane instanton wrapping the same cycle as the QCD 
stack. Its magnitude is set by ~ e~ 2nT and we show in table [1] the approximate 
magnitude of single- and double-instanton effects. In addition to the QCD cycle, 
there may be other cycles which instantons may wrap. We do not include these 
for two reasons: first, whether such instantons would generate a potential for 
the QCD axion is model-dependent^, and secondly, we can always arrange the 
model such that the QCD cycle is smallest and hence dominates the instanton 
expansion. We observe that the magnitude of the required cycle volume, and 
thus the magnitude of potential instanton effects, varies significantly with the 
string scale. If present, such D-instantons would generate a potential for the 
QCD axion. To compare their magnitude to that of QCD effects, we need an 
estimate of their contribution to the scalar potential. In this context we only 
care about terms containing a phase and so contributing to the axion potential. 
To this end, the relevant term from the scalar potential (123]) is 

V axion = e K (K? {diW(djK)W + c.c.)) . (62) 

9 It seems odd that such instantons could affect the QCD axion at all. However, if QCD is 
on cycle 3, and the axion 62 + 63 is fixed by effects on cycle 2+3, an instanton solely on cycle 
2 effectively generates a potential for the QCD axion. 
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Table 2: Magnitude of Axion Potentials from Superpotential Instanton Effects 



Euv 


10 8 GeV 


10 10 GeV 


10 12 GeV 


10 14 GeV 


10 16 GeV 


V 


10 20 


10 16 


10 12 


10 8 


10 4 


^1-instanton 


KT 79 M£ 


1(T 81 M 4 


10" 83 M| 


10" 85 M| 


10" 87 M 4 


^2-instanton 


10- 122 M| 


1(T 131 M 4 


10" 139 M| 


i(T 151 M 4 


1(T 155 M 4 



A superpotential instanton contribution e 27mT * generates a term 

V axton = "^^ e- 2 — 'cos(^). (63) 

The absolute magnitude of (163]) depends on the value of n, the internal volume V 
and the tree- level superpotential Wo- As we are looking towards phenomenology 
we also take m 3 / 2 = ^ ~ ITeV ~ 10 _15 Mp, as appropriate for gravity-mediated 
TeV-scale soft terms. Note that for models built around the KKLT scenario, we 
always have m s > Mqut and only the largest value of Eyy is achievable. In table 
[2] we give the internal volumes required for each UV scale, as well as the resulting 
absolute magnitude of l-,2- and 3-instanton superpotential corrections to the 
scalar potential. For the same reasons as above, we only consider instantons 
wrapping the QCD cycle. 

As well as superpotential effects, there are also nonperturbative corrections 
to the Kahler potential (the perturbative corrections to K, do not have an axionic 
dependence). While smaller, these are easier to generate - the instantons can 
have four fermionic zero modes rather than only two. A correction 

K = -2 ln(V) -»• K = -2 ln(V + e" 2 ™ T ) (64) 

will generate effects in the scalar potential at order 

W 2 

v SK ~ 9- e — 

and thus generate a potential for the QCD axion 6 of the form 

V&k cos{6 + a). 

Again assuming TeV-scale (visible) SUSY breaking, ^ ~ 1(T 15 , the magnitudes 
of such effects are shown in table [31 

The axion potential originating from QCD effects and relevant to the strong 
CP problem is 

Vqcd ~ Aq CD (1-cos(0)), 

with A QCD ~ 2 x 10- l9 M P and A 4 QCD ~ 10" 75 M 4 . We require QCD effects to be 
sufficiently dominant to be consistent with the failure to observe CP violation in 
strong interactions. Suppose we have a potential 

V = A(l-cos(#)) + ecos(# + 7). (65) 
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Table 3: Magnitude of Axion Potentials from Kahler Potential Instanton Effects 



Euv 


10 8 GeV 


10 10 GeV 


10 12 GeV 


10 14 GeV 


10 16 GeV 


V 


10 20 


10 16 


10 12 


10 8 


10 4 


^1-instanton 


1(T 99 M 4 


l(T 97 Mj, 


10" 95 M| 


1(T 93 M| 


1(T 91 M| 


^2-instanton 


10- 143 M| 


1(T 147 M 4 


1(T 151 M 4 


1(T 155 M 4 


1(T 159 M 4 



If A > e, the minimum is displaced from 9 = by 89 ~ 4. Observationally, 
\9\ < 10~ 10 , and thus non-QCD contributions must have absolute magnitude 
smaller than l(r 85 Mj,. By comparison with tables [2] and [3] it follows that in order 
for QCD instantons to dominate the axion potential, the instanton corrections 
to the Kahler potential are not constrained, while the 1-instanton superpotential 
correction may or may not be present depending on factors of 2ir and the precise 
value of the string scale. Higher instanton corrections to the superpotential are 
not constrained. 

To contribute to a superpotential (Kahler potential) an instanton must have 
at most 2 (4) fermionic zero modes, to generate J d 4 xd 2 9 and J d A xd 2 9d 2 9 terms 
respectively. In the absence of flux, there is a necessary condition on a divisor to 
generate a superpotential [21]: it must have holomorphic Euler characteristic one, 
Xg{D) = 1. In the presence of flux, this condition may be relaxed. The number 
of zero modes on an instanton, and thus its ability to appear in either the Kahler 
or superpotential, may also be affected by the presence of the stack of QCD 
branes wrapping the would-be instanton cycle. We shall not attempt to analyse 
this question for specific 'real' models, but by fiat simply assume the necessary 
instantons to be absent from the potential. In this regard it is encouraging that 
the number of instantons required to be suppressed is quite limited. 

4 The Axion Decay Constant 
4.1 Magnitude 

The last section was devoted to solving the strong CP problem: keeping an axion 
light while stabilising the moduli. However, even achieving this does not resolve 
all phenomenological problems. Given that a QCD axion exists, as indicated in 
the introduction there are strong bounds on the axion decay constant f a of equa- 
tion Q: 10 9 GeV < f a < 10 12 GeV. The lower bound, from supernova cooling, is 
hard. While the upper bound may be relaxed by considering non-standard cos- 
mologies, here we shall also treat this as hard. We want to estimate f a in some 
moduli stabilisation scenarios. In IIB compactifications, the axionic coupling to 
QCD arises from the clean and model-independent Chern-Simons coupling. How- 
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ever, to obtain the physical value of f a the axion must be canonically normalised. 
This depends on the Kahler metric, and in particular on where the moduli are 
stabilised. 

We assume we can write the Kahler potential as 

IC = JC(T l + T l ), (66) 

with /C real. This is true in perturbation theory, owing to the axionic shift 
symmetry, and any nonperturbative violations are small enough to be irrelevant 
for this purpose. In this case the kinetic terms for the axionic and size moduli 
do not mix. Noting that = JCjj, we have for any i and j 

= ^(28^8^ + 28^8%), (67) 

and the two sets of terms decouple. 

Let us first show that if both the overall volume and the cycle volumes are 
comparable to the string scale, then as expected f a > 10 16 GeV. Suppose an 
axion q is to be the QCD axion. The Lagrangian for this axion is 

Kad^dPa + ^jF a AF a . (68) 

For simplicity we have not included mixing terms: these will not greatly affect 
the discussion. 

The simplest toy model is that of a factorisable toroidal orientifold, with 
Kahler potential 

K = -ln((T 1 + f 1 )(T 2 + f 2 )(T 3 + f 3 )) 

= -ln(Ti + 7\)-lii(T 2 + r 2 )-lii(T3 + r3). (69) 

and Kahler metric 

/(Tx + fi)- 2 \ 

IC i3 = (T 2 + f 2 )- 2 . (70) 

V o o (r 3 + r 3 )- 2 / 

If we denote the axions by ci, c 2 and c 3 , the axion kinetic terms are 

-^Cid^Ci + ^2^c 2 ^c 2 + ^2^c 3 ^c 3 . (71) 

For definiteness, let us assume QCD is realised on cycle 1. There is no inter-axion 
mixing and the relevant axion Lagrangian is 

^d^Pd + ^ J F« A (72) 
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If we canonically normalise c[ = -^r, this becomes 

-d^d^ + ^^J F a A F a . 
In units where Mp = 1, the axion decay constant is 

f 1 



(73) 



If QCD is to be realised on this cycle, we need n ~ 12, and thus f a ~ 10 16 GeV. 
Going beyond this toy example, we recall that in general the Kahler metric was 
given by 



jc - = 



G 



3 ( k^jfot 3 ki mn t t kjpqt^t^ 



V 



V 2 



(74) 



If all cycles are string scale in magnitude, then ~ 0(1) and it is impossible 
to lower the axion decay constant substantially through canonical normalisation. 
The same conclusion applies: f a > 10 16 GeV. This conclusion is unsuprising: 
the axionic coupling to matter is a stringy coupling, and so we expect f a to be 
comparable to the string scale. If the string and Planck scales are identical, f a 
cannot lie within the allowed window. 

If we lower the string scale, phenomenological values for f a can be achieved. 
To analyse this, let us return to the toy model of fl52l) . We recall the Kahler 
potential was 



K 



The Kahler metric for this model is 



2 In ( (71 + fx) 5 - (T 2 + T 2 ) a - (T 3 + f 3 



(75) 
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\ 
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_3 1 9r 3 



(76) 



The axion kinetic terms are K^d^Cid^Cj 



At small volumes there is substantial 
mixing between the axions c\, c 2 and C3. However, in the limit V — > 00 with 
Ti ^ T2, T3, the Kahler metric takes the schematic form 



y-4/3 y-5/3 y-5/3 

y-5/3 y-1 y-2 
y-5/3 y-2 y-1 



(77) 



and is to a good approximation diagonal. The requirement t% 1 implies that 
QCD cannot be realised on branes wrapping cycle 1, as the resulting field theory 
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is far too weakly coupled. However, if T2 ~ T3 ~ 10 we may realise QCD by 
wrapping branes on one of these cycles. The resulting axion decay constant is 



a 




M P 



(9(1) 



M P . 



(78) 



Thus if V ~ 10 14 and r 3 ~ 10, the QCD gauge coupling is correct and the axion 
decay constant f a ~ 10 10 GeV lies within the narrow phenomenological window. 
Up to 0(1) factors, the string and Planck scales are related by 



Such a large volume corresponds to lowering the string scale to m s ~ 10 11 GeV. 
The lowered axion decay constant is easy to understand physically. f a measures 
the axion-matter coupling, which is an effect localised around the small QCD 
cycle. Thus the only scale it is sensitive to is the string scale, and so up to 
numerical factors f a ~ m s . 

The above is a very particular limit of moduli space, with one cycle taken 
extremely large while all others are only marginally larger than the string scale. 
It would thus be essentially a curiosity if it were not also the exact regime in which 
the moduli are stabilised in the compactifications of [TOl [TTj reviewed above. As 
the stabilised volume is exponentially sensitive to the stabilised dilaton, a priori 
the string scale can lie anywhere between the Planck and TeV scales. There is 
no difficulty, and no fine-tuning, in stabilising the volume so as to achieve an 
intermediate string scale. 

The above result on /„ is independent of whether an axion remains massless 
or not. As described in section I3T21 the simplest version of the scenarios of [T0l[TT] 
makes all axions far too heavy to solve the strong CP problem. In sections 13.41 
and 13.51 we have described the necessary modifications to this scenario such that 
a massless QCD axion will survive to solve the strong CP problem. Combining 
this with the above, we have for the first time given a procedure to stabilise all 
moduli while ensuring a QCD axion exists with a phenomenologically allowed 
value for f a . 

In itself this is interesting, as axions within the phenomenological window have 
always been hard to achieve in string compactifications. However, this scenario 
compels a further very interesting relationship between the axion decay constant 
and the (visible) supersymmetry breaking scale. 

4.2 Relation to Supersymmetry Breaking Scale 

We argued earlier that if a QCD axion is to be present in IIB flux compactifi- 
cations, QCD must be realised on a stack of D7-branes. We have also described 



g s Mp 



(79) 
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how to stabilise the moduli such that a QCD axion can exist with a phenomeno- 
logically allowed decay constant. 

By definition, any scenario of moduli stabilisation determines the moduli vevs 
and masses. However, in general much more information can be extracted. If 
the moduli potential breaks supersymmetry, it will generate soft supersymmetry 
breaking terms in the visible sector. In gravity-mediated scenarios, it is these 
that dominate the soft MSSM Lagrangian. It is an old and important problem 
to go from the moduli potential to the soft terms. The analysis of this problem 
was initiated in the heterotic string. As full moduli stabilisation in that context 
is hard, the relevant moduli potential was unknown. Progress was nonetheless 
achieved by parametrising supersymmetry breaking as S,T or U-dominated, de- 
pending on the moduli multiplet (dilaton, Kahler or complex structure) in which 
the dominant F-term occurred. This allowed a classification of supersymmetry 
breaking possibilities even in the absence of a full moduli potential. 

This question has now been revisited in IIB flux compactifications, where 
the moduli potential is under better control and full stabilisation is achievable 
(at least at the level of effective field theory). Instead of having to assume the 
structure of the F-terms, they may now be computed directly from the moduli 
potential. However, the problem is still subtle as the vacuum is originally AdS 
and the form of the (necessary) uplift to Minkowksi space can introduce extra 
contributions to the F-terms. In the KKLT scenario, the AdS minimum is su- 
persymmetric and all details of the supersymmetry breaking therefore lie in the 
uplift. Unfortunately this is the part of the construction under least control. Us- 
ing some particular assumptions about the uplift potential, (gravity-mediated) 
supersymmetry breaking in this scenario has been studied in |3Uj . 

For the exponentially large volume compactifications, the structure of super- 
symmetry breaking has been studied in (TT, 29J. The original AdS minimum 
is non-supersymmetric, and the structure of the F-terms is essentially inherited 
from the no-scale potential of [21J. While the uplift is again not well controlled, 
fortunately its effects are subdominant. This arises because the magnitude of 
the vacuum energy in the AdS minimum is O (t^), whereas from the F-terms it 
'ought' to be O (t^), the difference of course being accounted for by the no-scale 
cancellation. This implies that 'extra' F-terms required to lift the minimum to 
Minkowski space will be hierarchically smaller than those already present in the 
AdS minimum. 

No-scale breaking corresponds to Kahler-domination, with the dominant F- 
term associated with the multiplet controlling the overall volume [TTJ . The struc- 
ture of the soft terms on D7-branes for the no-scale models of [21] has been well- 
studied. Supersymmetry breaking is transmitted to the observable sector through 
gravitational interactions, and the 3-form fluxes present induce soft terms on the 
worldvolume theory of the D7 branes[36l 137] . Here we are only concerned with 
the results rather than the calculational details: the scalar and gaugino masses 
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are found to be 

m D7 ~ M D7 ~ m 3/2 = e ^ W „ f^. (80) 

If the internal volume V is exponentially large, it is the prime determinant of the 
scale of the soft terms. Other, more model-dependent, factors will also be present, 
but these will only be 0(1) effects relevant for the detailed structure but not the 
overall scale. As the volume also sets the scale of the axion decay constant, this 
implies a numerical relationship between these two quantities. From the fact that 

, M P , M P 

fa ~ and m soft 



it follows that up to numerical factors 

fa = ^M P m soft . (81) 

Thus in such models the axion decay constant is compelled to be the geometric 
mean of the Planck scale and the (visible) supersymmetry breaking scale. This 
is a striking result, as the two pieces of physics are a priori entirely unrelated. 
The origin of this can be traced to the intermediate string scale, which has been 
argued to have interesting phenomenological properties! 



5 Conclusions 

The purpose of this paper has been to investigate the conditions under which 
a QCD axion, ideally with a phenomenological value for /„, may coexist with 
stabilised moduli in string compactifications. This divides into two questions: 
first, how to stabilise the moduli such that a massless axion survives, and secondly, 
how to obtain allowed values of f a , 10 9 GeV < f a < 10 12 GeV. 

In the context of the first question, we have shown that the simplest version 
of many moduli stabilisation scenarios do not contain any light axions. We also 
have a negative result, in that supersymmetric moduli stabilisation is disfavoured: 
there exist no supersymmetric minima of the F-term potential with flat axionic 
directions. Even if AdS stability is present due to the Breitenlohner-Freedman 
bound, the tachyons must be removed by the time we are in Minkowski space. 
Performing this step requires a much greater technical understanding of uplifting 
AdS vacua to Minkowksi space than is currently available, and so it is unclear 
how relevant the original supersymmetric AdS solution is. 

This result is pure M = 1 supergravity and so makes no assumptions about 
the particular string model considered. It thus applies to all string compactifica- 
tions, and in particular shows that in many of the supersymmetric IIA flux com- 
pactifications considered recently the complex structure moduli sector is heavily 
tachyonic. 
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There is no no-go theorem on nonsupersymmetric minima with massless ax- 
ions, and so these may be preferred. In the context of the large volume compact- 
ifications of [101 [TTj . we outlined how to stabilise moduli while keeping axions 
massless. Here we had to rely on Kahler corrections that will become impor- 
tant as a cycle collapses to zero size. While unfortunately not much is known 
about these, our main requirement was simply that they exist. Clearly progress 
in determining the form of such corrections would be very interesting. We also 
specified the extent to which subleading higher-order instantons must be absent 
in order for a leading-order axion to solve the strong CP problem. 

We note this result also favours gravity-mediated supersymmetry breaking. If 
the moduli potential must break supersymmetry in order to solve the strong CP 
problem, then this suggests that supersymmetry should be broken at the string 
scale. Gravity mediation therefore always contributes to the visible soft terms 
and, unless the string scale is lowered to O(lOOOTeV), will dominate over gauge 
mediated effects. An intermediate string scale may then be preferred in order to 
obtain TeV-scale soft terms. 

In the context of the second question, the fact that f a is hierarchically lower 
than the Planck scale implies that compactifications with m s ~ Mp are unlikely 
to give allowed values for f a . In the compactifications of [10l [H] the string scale 
is hierarchically lower than the Planck scale. In these models, f a ~ M s and 
Msusy ~ An intermediate string scale therefore gives both 10 9 GeV < f a < 
10 12 GeV and visible susy breaking at O(lTeV). It is hard to find models with 
phenomenological values for f a , and so it is very interesting that in the above 
model this also implies TeV-scale supersymmetry breaking. We also emphasise 
that as very large volumes arise naturally in this model there is no fine-tuning 
problem in having msusY ~ O(lTeV). 

Moduli stabilisation and the landscape have received much discussion recently. 
Technically, the landscape large numbers of (9(10 500 ) arise from the very many 
ways of stabilising moduli. It is clearly necessary to find a handle for dealing 
with such numbers. One developed approach is statistical ([39] and references 
thereto). Without explicitly constructing examples of vacua, this aims at framing 
and answering questions about what is and what is not possible. However, it may 
be very difficult to identify the right vacuum: it has been recently argued that 
the problem of finding a vacuum with the right cosmological constant is NP hard 

A more general point argued here is that in the context of the landscape the 
strong CP problem may serve as an experimentum crucis. Assuming that the so- 
lution to the strong CP problem is a Peccei-Quinn axion and that string theory 
is a correct description of nature, this is a solution that is extremely sensitive to 
the physics of moduli stabilisation. Requiring an axion to remain (essentially) 
massless while all other moduli are stabilised is a technically clean problem di- 
rectly addressing the issue of vacuum selection. Indeed, as seen above imposing 
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this requirement directly rules out many scenarios of moduli stabilisation. The 
further condition 10 9 GeV < f a < 10 12 GeV is even more constraining: we have 
described one approach to this above and it would be very interesting to analyse 
others. 
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A Constraining Corrections to the Kahler Po- 
tential 

Here I discuss how to constrain the form of corrections to the Kahler potential. 
The tree-level Kahler potential used for our toy model was 



This geometry has one overall Kahler mode (Ti) and two blow-ups (T 2 and T 3 ). 
We specialise to our limit of interest T\ ^> r 2 , r 3 > 1, with t { = Re(Tj), where the 
resulting Kahler metric has the form (neglecting terms subleading in V) 



Any correction to the Kahler potential will also generate corrections to the Kahler 
metric fl83l) . As such corrections are perturbative, they may arise either from a' 
effects or loop effects. On physical grounds, we expect that such corrections will 
be subdominant to the tree-level metric in the regime - large overall volume, 
large cycle volumes and weak coupling - where both worldsheet and quantum 
corrections ought to be least important. 

We can apply this condition to restrict the form of potential corrections to K,. 
For example, consider the possible correction 



K = -2 In ((Tx + Tip - (T 2 + f 2 f 2 - (T 3 + f 3 ) 3/2 ) . 



(82) 




(83) 




(84) 
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with < a < 1. The 22 component of the corrected Kahler metric is 



/C 2 2 + SICoo = — ; — — o/n — • (85) 

2v^V 8V2r 2 3/2 V« 

As a < 1, the correction to the kinetic term would always dominate the tree-level 
term in the limit V ^> 1. This seems implausible as a large- volume limit ought 
to make the correction less, rather than more, important. 
A similar comment applies to a correction 

er 2 

/C + 5/C = -2m(V) + ^-, (86) 

which leads to 

/C 22 + <5/C 22 = = + (87) 

In this case, as we take r 2 large, the correction becomes increasingly dominant 
over the tree-level term. Given that large r 2 reduces both the curvature of this 
cycle and the gauge coupling on any brane wrapping it, we would again expect 
exactly opposite behaviour to occur. 

Finally, we could also consider the correction 

/C + 5/C = -21n(V) + ^, (88) 

with < a < 1. In this case 5K 2 2 is subleading to /C 22 in the classical limit. 
However, if we consider the 12 component we now have 

^ + ^ = -§# + Iy^ < 89 > 

and the correction again dominates in the large- volume limit. 
The correction used in the body of the paper, 

K + 5K = -2 ln(V) + ^ + (90) 

does not suffer from the above problems. As a correction to the Kahler metric, 
it gives 

/C 22 + 5/C 22 = * (91) 

2V2^V 16r 2 3/ V 

/C33 + 5/C33 = — ^=rr (92) 

Such corrections are suppressed compared to the tree-level term by a factor r _1 , 
i.e. g 2 of the field theory on the brane. Unlike those considered above, these 
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corrections are well-behaved (i.e. subdominant) in the classical limit, and there 
does not exist a 'bad' scaling limit. 

We note that for the case where Kahler corrections have been computed [41] . 
the correction does fall into this form. We only focus on the Kahler moduli 
dependence: the full expressions can be found in |JT]. The correction gives 

3 

/C + 8K = - ln(Ti + fx) - ln(T 2 + T 2 ) - ln(T 2 + T 3 ) + £ (93) 
with for example 

£n + ^ii = ^ + ^- ( 94 ) 

The loop-corrected Kahler metric is suppressed by a factor r-f = g 2 . 

The point of this discussion is that the form of possible corrections to the 
Kahler potential can be heavily constrained by the reasonable requirement that 
in a classical (large volume, weak coupling) limit, corrections to the metric become 
increasingly subdominant to tree-level terms: simply because Kahler corrections 
are very hard to calculate does not make us entirely ignorant of their form. 
In particular, denoting the 'small', blow-up moduli by Tj, these considerations 
exclude corrections of the form 

/C + 5K = -21n(V) + ^, (95) 

with a < 1, as in a classical, large volume limit there will be metric components 
whose correction dominates the tree-level term. 

This motivates our use of the correction (1901) : it comes from a reasonable 
assumption about the form of the corrections to the Kahler metric, and is con- 
sistently subdominant in the classical limit. 
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